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We propose a novel lattice calculation of spontaneous chiral symmetry breaking in QED model 
with 2+1 dimensional fermion brane. Considering the relativistic action with gauge symmetry 
we rigorously carry out path integral in Monte-Carlo simulation with Fermi-velocity relevant to 
effective coupling constant. We numerically show the evidence of spontaneous chiral symme- 
try breaking in strong coupling region with chiral condensate, low-lying mode distribution and 
Nambu-Goldstone boson spectrum in bare Fermi-velocty v = 0.1. This is a feasible study to 
investigate the phase structure of Graphene. 
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1. Motivation and background 

We consider QED model with 2+1 dimensional fermion brane (DQED3) |J^, ||]; in this model 
fermion has smaller velocity than speed-of-light c whose magnitude has been roughly estimated as 
c/300 Tight-binding approximation of Graphene []3|] (or similar material formed as honeycomb 
seat of atoms) has suggested that will possess a particular property which is described by a massless 
Dirac particle (quasiparticle) [Q, ||]. Although there are many analytic studies [7[ ||, ^] and 
numerical works jlO| , [Tl| ], it has not been successful to clearly describe the electric property of 
Graphene. 

Many models argue that there should be a signature of spontaneous gap generation due to 
strong dynamics which indicates the occurrence of semiconductor-insulator transition for sus- 
pended monolayer Graphene. However in the experiment there is no evidence of the gap of band 
structure between electron-hole, and thus it turns out to be permanently conductor. Such incon- 
sistency may be due to the suppression effect of silicon bases, whose dielectric permittivity £ 
appears in the denominator of effective coupling, a e (e) = 2a e / (1 + £)<1, and thus the impurity 
of Graphene might be cause of no phase transition. To verify this argument the identification of 
critical point of effective coupling constant from theoretical study is needed. Our motivation is 
clarification of phase structure of bQED3 model to understand the transition of Graphene and other 
kinds of honeycomb crystal. 



2. Brane QED 3 model 

As pointed out in [Q] the relativistic bQED3 model enables us to rigorously deal with path 
integral under gauge invariant formalism, 



S b QED 3 = f/ dtdx 3 (vE 2 + v l B 2 ^ + J dtdx 2 \fr iD t y t + iv(D x y x + D y y y 



¥ (2-1) 



with bare Fermi-velocity parameter. In this action we consider not only the Coulomb interactions 
but also induced magnetic interactions which are v 2 times weaker than electric field. Since we 
attempt to treat bQEDs model as local field theory (super-renormalizable), the theoretical param- 
eters, v and /3 , are affected by renormalization. The argument of perturbation theory in bQED3 



model is that the Fermi-velocity has a logarithmic divergence. The experiment [ |14] ] supports the 
renormalization effect of fermi-velocity of Graphene as suggested in perturbative bQED3 model 
[12, |l3|] rather than non-relativistic one [||, [7|, [g, |9[ [R], 11]. The effective coupling constant is 



defined as the modified form as a e = 1/ (471/3 v) and 2-loop analysis in approximated perturbation 
explicitly shows the UV fixed point in a e S> 1 []13|]. The check of existence of UV fixed point with 
non-perturbative method is interested in this work. 

The gap generation of Graphene is considered to be related with spontaneous chiral symmetry 
breaking (#SB) in bQED3 model. This is analogous to the second order chiral phase transition 
in Nf = 2 massless QED3 involving the mass gap [ PI 16]. In Eq.([Q|) quasiparticle field y has 



4 spinor component, \j/ = (Ya + i¥o + >Vct _ >Vv _ ) which corresponds to three kinds of Graphene 
symmetry; Dirac valleys (degenerating ground energy) (±), sublattice symmetry (A,B) and spin 
of carbon atoms (a). The gamma matrix is given as the tensor structure; y t = Oq ®hxi, Ji = 
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— i<72 <g> (T;(i = x,y,z) in which the first part is degree of freedom (DOF) in valley times sublattice 
and the second one is DOF in spin rotation. Regarding the Graphene symmetry as "flavor" U(4) 
symmetry in DQED3, whose 16 generators are represented as {1,75,^, [73, 7s] /2} <S> Oj =t ^^ z with 
75 = iYtYxYyYz* we define 75 = diag (1,1, — 1, — 1) "chiral" projection corresponding to left and right 
chirality (valley chirality) as Ul(2) x Ur(2). The chiral condensate (yy) is use d to be the order 
parameter of 75 #SB as well as QCD. 



3. Lattice calculation of brane QED3 model 



In order to carry out ab initio calculation of DQED3 model, we implement the staggered-type 
fermion action including the non-compact U(l) gauge action; 



I 

n=(x,y,z,t) 

I 

m=(x,y,t) 



j3v£(V 4 e ; -(«) - V;0 4 («)) + /3V 1 £(Vj0 y (n) - VMn)) 

i i,j 

£ Hi ( m )£ ( m ) {Ui(m)x(m+T) - l/J (m)%(m -i)}+M% {m)% (m) 



(3.1) 



(3.2) 



with differential V,(m,«) = 5 m ,, 1+ i — 5 TO>n , link variable E7,-(m) = exp(i0(m)) and Kawamoto-Smit 
phase factor TJ,-(m) = IIfc=i ( — !)*• Conveniently the staggered action in 2+1 dimension has U (2) x 
U (2) (flavor-) chiral symmetry in the continuum limit, and hence we can perform Hybrid-Monte- 
Carlo (HMC) simulation straightforwardly [jlO[ [H], §]. 

Here we show the numerical results of chiral condensate, low-lying mode distribution and 
hadronic spectrum in N s x N t = 40 2 x 20 and N z = 8 lattice with fixed v = 0. 1 as a function of 
fermion mass M and coupling constant j8. For U(l) field we use periodic boundary for spatial 
and temporal directions, and for fermion field is periodic boundary for spatial direction and anti- 
periodic boundary for temporal one. This simulation realizes the lower temperature system thanks 
to a rescaled temporal size by velocity. To avoid the autocorrelation we use every 20 HMC tra- 
jectory per configuration, and the error analysis adopts the Jackknife(JK) method with 10 bin size. 
Chiral condensate is estimated by 100 noise sources in each configurations. 



4. Chiral symmetry breaking in brane QED3 model 

4.1 Chiral condensate 

We first show the chiral condensate and chiral susceptibility; 

a^{ X x) = (l^D-\ mi m)), Zm ^|^ = /(£D- 1 (m,m)) 2 -l£|D- 1 (m,0)| 2 ) (4.1) 

m m m 

in several j8 and M, where D 1 denotes the inverse of staggered fermion matrix. Figure [I] shows 
a and Xm nave clear dependence of both bare effective coupling a = 1/(4jtj3v) and mass M. a 
drastically grows up at j8v ~ 0.05 - 0.06, and decreasing M this growth is further developing. It 
seems that the critical point exists around j8v ~ 0.05 - 0.06 which corresponds to a c ~ 1 .3 - 1 .6, and 
in a < a c G is close to zero while in a > a c G remains in finite value. The chiral susceptibility is 
clearly shown to be significant M and l/(j8v) dependence similar to critical transition. In a c ~ 1.3 
-1.6 there will be a singular point for Xm in M = limit, which indicates that a c is expected to be 
discontinuous at the critical point. 
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Figure 1: The dependence of the inverse of effective coupling constant 1 / Q3v) for chiral condensate cr (left) 
and chiral susceptibility % m (right). Different symbols denote the points in the different fermion masses. 



4.2 Low-lying mode distribution 

Investigation of low-lying mode of massless Dirac operator (Dirac kernel) is helpful to quantify 
renormalized chiral condensate £ through the Banks-Casher relation; 

I/7r = lim hmp( £ ), p(A) = U^8(X - X n )\, (4.2) 

in chiral broken phase. Spectral density p (X ) is as a function of the eigenvalue X n of Dirac kernel 
(here we consider Dirac kernel in finite volume V, and thus distribution of A is a discretized distri- 
bution.). The Banks-Casher relation expects that spectral density shows a constant distribution near 
zero eigenvalue. In finite volume, according to the discussion of random matrix theory (RMT), the 
"hard edge" [ |l7| ] of p appears at X = with width 1 / (VL). On the other hand in symmetric phase 
the spectral density p(X) ~ X d ~ x (d denotes dimension of fermion), at weak coupling regime. This 
different behavior is used to not only distinguish the phase in effective coupling constant as well as 
Figure [T] but also quantitatively estimate the £ in constant region of p near X ~ 0. In the calculation 
of low-lying eigenvalue we evaluate 75 different eigenpairs (staggered-type Dirac kernel has a pair 
of eigenvalue with different sign). 

In Figure ^| we can clearly see a significant change of the low-lying distribution of p above 
and below a c which is observed in Figure ||. Below a c (j8 ~0.4), there is clear plateau starting 
from X = 0.002 in which spectral density is consistent with one in different volume. Near zero 
point there is also expected "hard edge" whose width becomes narrow increasing volume as in the 
chiral broken phase. When decreasing the a, this "hard edge" disappears and in a < a c p becomes 
monotonically growing. Fitting in 0.002 < X < 0.007 with constant function, we obtain £ in each 
sea fermion mass. 

4.3 Nambu-Goldstone boson spectrum 

In order to confirm £ obtained in spectral density, we attempt to reproduce this from measure- 
ment of the spectrum of Nambu-Goldstone (NG) boson particle in bQED3. In NG theorem, when 
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Figure 2: Left panel shows the spectral density at different j3 which are above, close and below a c . We also 
compare different spatial volume. Right panel is a zoomed-up right panel near X = 0. This is result in sea 
fermion mass M sea = 0.001. 



spontaneous #SB occurs, NG bosons appear as an asymptotic one-particle state coupled with axial- 
vector current (or pseudoscalar via Ward-Takahashi identity) operator. If a > a c is in the chiral 
broken phase, four NG bosons exist as well as pion (and eta) meson in QCD (note that in bQED3 
model, since there is no anomaly, flavor singlet particle can be regarded as massless NG boson.) 
If this picture is true, we can extract the NG boson amplitude and mass from pseudoscalar (PS) 
correlator in large distance separation between source and sink point. After taking large separation 
of x direction we expect that PS correlator which is defined as the exact NG boson operator in 
staggered-type fermion approaches to the following form: 

limGpsW = \Z NG \ 2 (e- m ^ x + e- m ^ N ^)/(2m NG ), / NG = 2MZ PS m NG /2 , (4.3) 

where /no is a corresponding quantity to pion decay constant in QCD although its mass dimension 
is 1/2 in 2+1 dimension. Here we set the source operator in the origin as point source. 

Figure || shows the clear shape of exponential function of Gps, and effective mass plot explic- 
itly illustrates that in a > a c there is plateau in x > 9, however in a < a c plateau is not observed. 



This result indicates that NG boson state following in Eq.([4.3[) appears as a consequence of spon- 
taneous #SB in a > (Zc- 



Using the value obtained from fitting Gps with function in Eq.(|4.3|), we show the lattice results 
of /ng and m NG /M in each NG boson mass squared in Figure f|. To evaluate the chiral condensate 
from these results we can use the GMOR relation; I = / m|/(2M) as in the case of QCD where fo 
is a value of /ng in M = 0. According to the loop calculation in chiral perturbation theory (ChPT) 
in 3-dimension, we set the fitting function 

/NG = /o(l+>W(47r/o 2 ))+«4 G , ml/M = 2B{l-m NG /(47ift))+dm^ G , (4.4) 

where fo, B, c and d are fitting variables. We introduce c and d as correction terms to higher order 
effect than LO ChPT. We attempt to compare the extrapolated result with different fitting range of 
m 2 K , and including the data at M = 0.005 to extrapolate M = limit we use a function in Eq.([k4]) 
besides using data up to M = 0.0025 we exclude c and d terms. In Figure |4] we see that in a > a c 
(/3 = 0.4, 0.45, 0.5) the chi-square fitting with both NLO ChPT and NLO CHPT + linear term 
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Figure 3: Gps(x) (left) and effective mass plot (right) at different j3 in M = O.OOl(circle), M 
0.0025(square) andM = 0.005(diamond). The straight lines are fitted result. 



Table 1: Fitting results of chiral condensate obtained by spectral function via Banks-Casher relation and 
ChPT with NLO and NLO+linear. We also describe the value of ^ 2 /dof. The error is only statistical one. 



p 


Z (spectral) 


£ (NLO ChPT) 


% 2 /dof 


£ (NLO ChPT+linear) 


% 2 /dof 


0.4 


0.2433(48) 


0.2587(8) 


12.2(4.1) 


0.2657(15) 


6.1(2.9) 


0.45 


0.2140(45) 


0.2243(7) 


4.1(2.3) 


0.2275(13) 


2.7(1.9) 


0.5 


0.1694(51) 


0.1807(6) 


1.5(1.4) 


0.1766(13) 


2.9(2.0) 



works better than in a < a c . This is also consistent with picture that spontaneous #SB above a c 
occurs in accompany with NG boson. In table [l] we show preliminary results of comparison of £ 
with three different ways in above a c . These are consistent results within l-2a under a few % 
accuracy. 



5. Summary and discussion 

We perform the ab initio calculation of QED model with 2+1 dimensional fermion brane 
(bQED3) in Monte-Carlo method. Due to preserving gauge symmetry we take into account not 
only coupling constant but also fermi-velocity as theoretical parameters. In this proceedings we 
numerically show the strong evidence of ^SB and critical coupling constant from three kinds of 
way; chiral condensate, spectral density and NG boson spectrum. Although we fix the "bare" 
velocity parameter in 0.1 which is relatively larger than naive estimate (v ~ (?(10~ 3 )), we have the 
consistent result with expected in #SB phenomena. This is feasible study to search the applicability 
to study of phase structure of Graphene. 

The calculations were performed by using the RIKEN Integrated Cluster of Clusters (RICC) 
facility. This work is supported by the Grant-in-Aid of the Japanese Ministry of Education (No. 
20105002, 23105714(MEXT KAKENHI grant)). 
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Figure 4: M dependence of /ng (left) and m\/M (right). Solid line denote the LO ChPT in the range of 
0.001 < m^Q < 0.0025 and dashed one denotes the fitting function in LO ChPT plus linear term in the range 
of 0.001 <m^ G < 0.005. 
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